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AUSM 系列 算法 对 比 研究 及 背景 太阳 风 初 步 应 用 
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摘 要 磁 流体 力学 数值 模拟 是 研究 日 地 物理 学 现象 的 一 个 重要 手段 . 对 比 三 种 AUSM 算法 , 即 AUSM, AUSM+ 
和 AUSMPW--, 结合 HDC 磁场 散 度 消去 方法 计算 多 维 MHD 问题 的 性 能 . 通过 分 析 三 种 算法 计算 Rotor 算 例 
和 Orszag-Tang vortex 算 例 的 结果 发 现 , AUSM+ 算法 的 性 能 最 好 . 进一步 使 用 AUSM-- 算法 基于 6 片 网 格 构 
造 模拟 了 日 园 结 构 , 计算 结果 表明 这 种 算法 能 够 正确 计算 出 日 早 的 大 尺度 结构 . 对 于 日 园 结 构 模 拟 块 ， HDC 方法 
能 够 较 好 地 控制 磁场 散 度 误差 . 
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Abstract Numerical simulation of MHD (Magnetohydrodynamics) is an important method to 
study solar-terrestrial physics phenomena. In this paper, three AUSM algorithms, i.e., AUSM, AUS- 
M+, AUSMPW+, were compared in conjunction with the HDC method to eliminate the divergence 
of the magnetic field. The AUSM- algorithm is found to be better than the AUSM algorithm and 
the AUSMPW+ algorithm through the results of Rotor example and Orszag- Tang vortex example. 
Further, the AUSM+ algorithm is used to simulate the coronal with a six-component grid system. 
'The results show that this algorithm can correctly calculate the large-scale structure of the corona. 
Also, the HDC method can maintain the divergence-free constraint on the magnetic field. 
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0 引言 


H3 zx HH PECES SI HA, TER, i 
层 、 电 离 层 等 物理 性 质 不 同 的 空间 区 域 ,， 磁 流体 力 
学 (MHD) 数值 模拟 是 研究 日 地 物理 学 现象 的 一 个 
重要 手段 , 被 广泛 应 用 于 日 冕 磁场 结构 1、 背景 太 
阳 风 结构 9. BV (CME) 事件 向 、 太 阳 
风 - 磁 层 相 互 作用 等 研究 Pl. 随 着 基于 磁 流体 力学 数 
值 模拟 的 空间 天 和 气 预 报 方法 投入 运行 回 , 空间 天 气 
预报 由 经 验 预报 向 数值 预报 发 展 . 

在 MHD 数值 模拟 过 程 中 , 一 般 需 进行 特别 处 
理 以 使 磁场 散 度 误差 (V . B = 0) 保持 为 零 I. H 
前 主要 的 磁场 散 度 误差 处 理 方法 有 4 种 . (1) 泊 松 校 
正法 9, 通过 求解 泊 松 方程 计算 磁场 中 与 散 度 误差 
相关 的 部 分 , 并 将 其 从 磁场 计算 结果 中 消去 ; (2) CT 
(Constrained Transport) 方法 °, 通过 求解 电场 旋 度 
形式 的 磁感应 方程 , 使 磁场 散 度 误差 不 随 计算 增加 ; 
(3) Powell 的 8 波 模型 DO, 在 守恒 形式 的 MHD 方程 
组 右 端 添加 正比 于 V-B 的 源 项 , 使 得 磁场 散 度 误差 
可 以 随 着 流体 运动 被 输 运 出 去 ， 避免 散 度 误差 累积 ; 
(4) Dedner 的 HDC (Hyperbolic Divergence Clean- 
ing) 方法 M, 通过 在 MHD 方程 组 中 添加 一 个 耦合 
磁感应 方程 和 V . = 0 的 方程 , 使 得 MHD 方程 组 
变 为 9 个 方程 , 在 不 破坏 MHD 方程 守恒 性 及 不 增加 
额外 计算 步骤 的 情况 下 消除 磁场 散 度 误差 . 

在 日 地 空间 物理 学 的 MHD 模拟 研究 中 大 都 采 
用 有 限 体积 类 数值 方法 . 求解 有 限 体积 方法 控制 
体 界面 通 量 时 ,迎风 格式 根据 方程 组 特征 波 的 传播 
方向 计算 通 量 , 能 够 正确 反映 物理 信息 的 传播 方向 ， 
获得 较 好 的 间断 捕捉 效果 . 迎风 格式 可 以 分 为 FDS 
(Flux Difference Splitting)!? 和 FVS (Flux Vector 
Splitting)?! 两 类 FDS 通过 界面 上 近似 黎 曼 问 
题 的 解 得 到 界面 上 的 通 量 ， 最 典型 的 FDS 类 格式 
是 Roe 格式 ; FVS 格式 则 直接 在 界面 上 分 裂 通 量 
的 Jacobi 矩阵 得 到 包含 迎风 信息 的 数值 通 量 , 典型 
的 FVS 类 格式 有 Van Leer?! 和 Steger-Warmingl 
格式 . 为 了 集合 FDS 离散 格式 的 高 精度 与 FVS 
离散 格式 的 稳定 性 US, 构造 了 混合 通 量 分 裂 格 式 ， 
具有 代表 性 的 有 AUSM (Advection Upstream S- 
plit Mach number), LDFSS (Low-diffusion Flux- 
splitting Scheme)!!3 和 CUSP (Convection Upstream 
Split Pressure)!!9] 等 . 
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AUSM 算法 是 迎风 格式 的 一 种 , 在 可 压缩 流体 
力学 (HD) 数值 模拟 中 取得 了 较 好 的 效果 . AUSM 
算法 的 核心 思想 是 把 方程 组 的 流通 量 分 为 流动 项 和 
压力 项 ， 然 后 根据 网 格 界面 的 马赫 数 计算 流通 量 . 
这 种 算法 不 需要 计算 通 量 Jacobi 矩阵 的 特征 值 和 
特征 向 量 , 方便 易 行 . 作为 从 计算 流体 引进 的 算法 ， 
1999 年 Agarwal 和 Augustinus? 首先 利用 原始 的 
ASUM 算法 计算 一 维 Brio-Wu 激 波 管 问题 虽然 
得 到 成 功 应 用 但 是 在 强 激 波 处 会 出 现 数值 振荡 ， 为 
了 提高 AUSM 算 法 的 稳定 性 ，2007 年 Pan?! 在 非 
结构 网 格 下 使 用 AUSMT+ 算法 计算 了 一 维 的 Brio- 
Wu 激 波 管 和 二 维 nozzle 流动 问题 , 取得 较 好 的 效 
AR. 为 使 AUSM 系列 算法 更 有 效 地 捕获 间 断 ，2009 
年 Hanba 将 气体 动力 学 改进 的 AUSMPW+ 算法 引 
入 磁 流 体力 学 求解 激 波 问 题 , 捕获 能 力 增 强 . 

本 文 介绍 了 AUSM, AUSM+ 和 AUSMPW+ 三 
种 AUSM 数值 格式 构造 和 控制 磁场 散 度 误差 的 HDC 
方法 , 为 了 对 比 这 三 种 算法 在 计算 多 维 MHD 问题 时 
的 性 能 , 计算 了 二 维 Rotor 和 Orszag-Tang vortex 两 
个 标准 算 例 , 并 使 用 计算 效果 最 好 的 AUSM+ 算法 
HAT HAA. 


1 数值 算法 
1.1 ”控制 方程 


采取 无 量 纲 化 变量 ， 控 制 方程 采用 的 是 如 下 和 守 
恒 形 式 的 单 流体 理想 MHD 方程 : 


à 
S V (ov) - 0, (1) 
pv 1 1 E 
sep ENS [wv + (ot)! =0, 
(2) 
2E +V [m z-Bh)v- 
=w BB = 0, (3) 
ƏB 
a tV (v:.B—-B.v)= (4) 
同时 , 磁场 遵循 高 斯 磁 定律 : 
V. B-—0. (5) 


其 中 , p 为 密度 , v 为 速度 , p 为 气体 压强 , B 为 磁 感 
应 强度 , uo 为 真空 磁 导 率 , y 为 气体 比 热 比 ,7 为 单 
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位 张 量 , E 为 能 量 , 有 


-P lu, lp 
uar iea” T : 
1.2 AUSM 系列 算法 
理想 MHD 方程 的 守恒 形式 可 以 写 为 


Qt OL (6) 
Ot Or Oy Oz ` 


AF, U 为 流动 量 , F 为 zx 方向 流通 量 , G y y 方向 
流通 量 , H 为 z 方向 流通 量 ， 


p 
pu 
pu 
w 
Bes t 
Bz 
By 
B; 
E 
pu 
pu? +p — Bz 
puv — Br By 
uw — B4B; 
F- » . 
uB; — uBzr 
(7) 
uBy — vB; 
uB, — WB 


(E + piu — Bs(uBz + vB, + wB;) 
半 离 散 化 的 守恒 型 三 维 磁 流体 方程 的 形式 为 
OU FT1127 — Fi—1/2,j,k 
72 rm * 
0y 
Hi jk+1/2 — Hijk-1/2 — 
eo n 0. (8) 


AP, Fii Fi-1/25j,k> Gage Gij-1/2,k> 

Hi jkay; Hij-ij2 为 控制 体 界面 上 的 数值 通 量 . 
在 计算 数值 通 量 时 ,AUSM 将 通 量 分 解 为 流动 

项 FC 和 压力 项 FP. 以 z 方 向 上 的 通 量 FP 为 例 ， 


F = FF 4 F? = uy + P. 
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其 中 ， 
eer Mu sal 

pu m - B; 
pu —B&y4 Dy 

rA pw pP- —B,DB; 
Ba -uB, |’ 
By —vB, 
B, —w By 

E +p —B,(v - B) 
l, p2 2 2 
p =p+3(B + B? + B?). (9) 


在 AUSM 系列 算法 中 , 不 同 算法 对 界面 上 的 通 
量 采 用 的 计算 方法 有 所 区 别 , 下 面 分 别 了 予以 介绍 . 
1.2.1 AUSM 算法 
Fijo, ausm = Ce (Mi Yr + Mg bg) 
(pi Pr + pg Pn). (10) 


式 中 , 1/2 为 网 格 界面 处 流通 量 , (L,R) 为 界面 处 重 
构 的 控制 体 变量 ， 


ML,R = Uu 
Cf 
1 
IM +1}, |M|<1, 
TA | ; 
了 c M|), |M] >al; 
ÍM +1)?@mM), |M|<1, 
z! + sgn(M)], |M| > 1; (11) 
c; = min (ct L, cg); (12) 


p peip a 
G&L = 4 z| H 
2 PL 


(um +Bz Duy B (ER a 
pL [7] 


- fp [et Pe Be 
CR 二 | 一 一 一 一 一 一 一 十 
2 pL 


(s t Bn: uy N er ia] ` 
PR PR | 
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By, = 51B 为 法 向 磁场 . 
1.2.2 AUSM+ 算法 
AUSMT+ 算法 与 AUSM 算法 的 区 别 在 于 构造 
, PE 时 引入 了 o 和 6 项 , 此 时 界面 的 流通 量 


M 

为 

Fi ja,AUSM+ = er(Mg |g-iys UL + Mg |g-iysUR)- 
(Pr |o-3/1e Pr + Pg |o-3/16 PR). 


(15) 
ME 和 PE 的 取 值 分 别 为 
EM 1 BI? -12 

Mis IM| & 1; (16) 

3(M +|M)), IM|>1. 

HU x 1? (2mM) € aM(M? — 1)? 
Pe IM| <1; 

ju :sgn(M), |M|» 1. (17) 


1.2.3 AUSMPW+ 算法 

AUSMPW+ 算法 在 AUSM+ 算法 的 基础 上 增 
加 了 压力 加 权 项 , 重新 构造 ME, 得 到 ME. 此 时 界 
面 流 通 量 


1 /2,AUSMPW+ = 
ce(Mi |p=1syr + Mg lg-iysR)- 
(Pf |a=3/16PL + Pa|a=3/16PR). (18) 
其 中 ， 
M: = Mj + Mg 2 0, 
w): (1+ fr)— fz], 


Mf = M$ Mj 
Mg = Mgw( + fn): 
M; = M} +M} « 0, 
Mf = Miw(- fr), 
Mg = Mg + M#[(L — v) ^ fr) — fn]. 
(19) 
AF, f,w 为 与 压力 相关 的 系数 , 即 


Pt, L, Pt, R 
一 一 一， pus 
ÍL,R = Pus 


0, pt,s = 0. 


Dus = Pi pL + PR PR. (20) 


Chin. J. Space Sci. 空间 科学 学 报 2015, 35(4) 
E PaL Pony? 
w= 1- min (Ese, PE) ; (21) 
Pt,R PuL 


1.8 ”磁场 散 度 误差 处 理 

为 了 限制 磁场 散 度 误差 的 增加 ， 使 用 Dedner 
的 HDC 方法 叫 处 理 磁场 散 度 . HDC 方法 通过 变 
E 2 将 磁场 散 度 方程 与 磁感应 方程 联系 起 来 , 即 


E+. B'—-B.v')+Vy=)0, 
D(g) 4 V. B «0. (22) 


AF D 为 线性 差分 算 符 , D 取 值 不 同 , 可 得 到 不 同 
的 修正 方法 P9741 


这 里 取 文献 [24] 中 的 
| 109 9 


Lo 


A (24) 为 变量 的 对 流 - 耗 散 方程 , 磁场 散 度 误 
差 随 变量 yp 被 输 运 到 计算 边界 , 在 输 运 过 程 中 被 耗 
散 减 小 . 式 (24) F, eu 为 输 运 速度 , cp 为 耗 散 率 . 为 
了 不 对 计算 的 时 间 步 增加 额外 的 限制 , 将 cn 取 为 
CFL 


Ch TUA max(1/d)' 


[du ean] 


这 里 CFr = max d 为 网 格 大 小 . 


按照 文献 23], 耗 散 率 
Cp AR 


此 时 , 控制 方程 式 (1)~(4) 变 为 
a OF 0G ôH 


E Ue que e een 
Do" inr 2 ^ ee 
pv 
p 
1 
jo pov - I(p- ;IBP) - B.B 
à 
ce eas v.B—- B.v4 Io 2 
E 1 
(E »* zIBIP)v - Bw: B) 
d 2 


2 
cL B 
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0 
0 
(25) 
0 
2 
二 
ce g 
界面 上 的 流通 量 相应 变 为 
ori 1 n n 
Fij2 = Fija + 3 (XB,L t XB,n) (26) 


AP, FU 为 前 文 推导 的 式 (10) AUSM 系列 算法 、 
式 (15) AUSM+ 或 式 (18) AUSMPW+ 的 界面 流通 
量 FE. 

VÀ x Jr Afi], 式 (26) 通 量 的 具体 表达 形式 为 


[e] 


P- 一 (Br T Bs, n)u/2 , 
一 (Br Te By,R)v/2 
一 (Br + Bz r)w/2 


一 (Br F Bz R)(v : B)/2 


XB-— (27) 
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1.4 EHHA 
重 构 方式 这 里 选择 了 MUSCL 方法 P3, Vp x jr 
向 为 例 , 重 构 过 程 可 以 表达 为 


1 
UL. —Ui jk t n s K) 


min mod (AUi. 172,5; bAU; 1/2, k)- 


1 : 
id + K) min mod (AU; 475,5 x, 0AU; 12,5); 
(28) 
1 
Uj S mU - d 7 K) 


2,9, 
min mod (AUi 3/2, bAU; 1 5k)— 


1 : 
n5 T K) min mod DAL x bAU; 1/2, k)- 
(29) 


式 中 ， 


AUDiHlya k = Ui+1,j,k — Uii, 
K = -1,0,3; b= 0.9. 


1.5 “时间差 分 
对 时 间 采 用 了 二 阶 Runge-Kutta 离散 法 . 设 半 
离散 化 方程 为 
E = Lr (u). 
式 中 , Lh(w) 为 空间 离散 算 子 . 
二 阶 精 度 时 间 差 商 格式 为 


ut — u? 4 AtLs(u?), 


1 1 
ut = (up u() 十 gAtL(uC)). — (30) 


j j 


2 数值 算 例 


为 了 验证 和 对 比 三 种 AUSM 数值 算法 的 性 
能 , 用 这 三 种 算法 分 别 计算 了 二 维 Rotor 算 例 和 二 
维 Orszag-Tang vortex 算 例 . 
2.1 二 维 Rotor 算 例 

Rotor 问题 26] 被 建议 作为 一 个 检验 Torsional 
Alfven Waves 运动 的 标准 算 例 . 这 个 问题 假设 开 
始 时 在 静止 流体 内 部 有 一 个 稠密 气体 形成 的 快速 旋 
转 圆 柱 体 . 本 文 取 二 维 笛 卡 儿 坐标 系 下 ，(2x,y) € 
[—0.5, 0.5]2. 

该 问题 的 初 值 为 


D=1，DB=5/V4r， B,= B= 0, 
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p—1l1-9f(r, r= +y. 
u = —2f(r)y/0.1, v =2f(r)x/0.1, 
w=0.0, r<0.l; 
u = —2f(r)y/r, v =2f(r)z/r, 
w=0.0, r>o0.l. 


1, 7 < 0.1; 


200 
f(r) = = (0118-7) 0.1 <r < 0.115; 


0, r>0.115. 


(a) AUSM 


0.2 04 
(b) AUSM+ 


-0.4 -02 0 02 04 
(c) AUSMPW+ 


图 1 AUSM 三 种 算法 对 Rotor 算 例 计算 得 到 的 
B. 等 高 图 (t = 0.15) 
Fig.1 The x direction magnetic field of three AUSM 
algorithms to solve the Rotor problem (t — 0.15) 
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绝热 系数 y = 1.4, 函数 f(7) 通过 让 2 个 量 在 
从 内 部 状态 (r < 0.1) 到 外 部 状态 (r > 0.115) 的 一 
个 小 距离 ôr = 0.015 内 , 以 线性 方式 平滑 圆柱 体 表 
面 密 度 和 速度 的 不 连续 性 . 网 格 数 为 200 x 200, 边界 
为 自由 输出 . 图 1 和 图 2 给 出 了 AUSM 三 种 算法 计 
算 Rotor 算 例 所 得 压强 和 B... 

在 初始 速度 作用 下 ,计算 区 域 中 的 流体 开始 旋 
转 , 同时 在 离心 力 的 作用 下 向 外 运动 , 使 得 中 心 区 域 
压强 下 降 . 对 比 图 1 和 图 2 所 示 的 三 种 算法 计算 结 


(a) AUSM p 


-0.4 -02 0 02 04 
(b) AUSM+ 


02 04 
(c) AUSMPW+ 


-0.4 -0.2 0 0.2 0.4 
x 


图 2 AUSM 三 种 算法 对 Rotor 算 例 计 算得 到 的 
p 等 高 图 (t — 0.15) 
Fig.2 Pressure field of three AUSM algorithms to solve 
the Rotor problem (t — 0.15) 
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果 可 知 , 三 种 算法 均 正 确 计算 出 Rotor 问题 的 结构 ， 
但 是 AUSM 算法 计算 结果 中 的 等 值 线 不 平滑 ， 
AUSMPW+ 算法 结果 中 的 等 值 线 出 现 比较 明显 的 震 
H, 而 AUSM+ 算法 的 计算 效果 较 好 . 
2.2 Orszag-Tang vortex 算 例 

由 于 Orszag-Tang vortex 问题 P7) 在 vortex 形 
成 时 有 多 种 激 波 的 相互 作用 , 常 被 用 来 作为 考核 磁 流 
体力 学 算法 有 效 性 的 标准 算 例 . 本 文 把 Orszag-Tang 
vortex 问题 的 初始 条 件 给 定 在 (v, y) € [0,27]? 的 一 
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图 3 AUSM 三 种 算法 对 Orszag-Tang vortex $345] 
计算 得 到 的 p 等 高 图 (t = 3) 
Fig.3 Pressure field of three AUSM algorithms to solve 
the Orszag-Tang vortex problem (t = 3) 
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个 正方 形 区 域内 , 有 


p(z, y,0)=7, 
v(z, y, 0) = sin z, 


p(z. y,0) = 6", 
u(x, Y, 0) = —sin y; 


B,(z,y,0)— —siny, By(x,y,0) = sin 2z. 


其 中 , 7 = 5/3, 采用 周期 性 边界 条 件 , 网 格 数 为 200 x 
200. 


图 3 和 图 4 给 出 了 三 种 算法 计算 Orszag-Tang 


(a) AUSM Bx 
6 
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(b) AUSM+ 
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(c) AUSMPW- 


图 4 AUSM 三 种 算法 对 Orszag-Tang vortex 算 例 
计算 得 到 的 B. 等 高 图 (t = 3) 
Fig.4 The x direction magnetic field of three AUSM 
algorithms to solve the Orszag-Tang 


vortex problem (t — 3) 
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vortex 算 例 的 压强 和 2z 方向 磁场 的 等 值 线 图 .对比 
其 他 文献 的 计算 结果 P7] 可知 , 三 种 算法 均 正确 计算 
出 Orszag-Tang vortex 算 例 的 总 体 结构 . AUSM+ 得 
到 的 效果 比较 好 , 而 其 他 两 种 算法 结果 的 等 值 线 都 出 
现 了 锯齿 状 非 物理 结构 . AUSM 的 耗 散 较 大 , 会 忽略 
一 些 细小 结构 . 图 5 为 三 种 算法 对 vortex 算 例 计算 
所 得 压强 . 从 图 5 可 以 看 出 , AUSMPW-- 的 捕获 能 
力 强 , 在 x = 2.6, 3.4 均 处 捕捉 到 激 波 , 但 是 同时 振 
荡 也 较 大 . 


3 AUSM+ SüETE H aaith 
数值 模拟 中 的 初步 应 用 


通过 前 文中 算 例 计算 结 果 的 对 比 , AUSM+ 算法 
能 够 得 到 比较 好 的 结果 , 因此 本 文 将 这 种 方法 进一步 
应 用 于 日 网 结构 的 数值 模拟 ， 计 算 区 域 选 取 为 从 太 
阳 表 面 到 16 Rs( 太 阳 半 径 ) HER, 采用 文献 [28] 的 
方法 , 将 球面 划分 为 相同 的 6 个 网 格 区 域 采用 并 行 
计算 ， 该 网 格 系统 有 效 解决 了 极 区 网 格 收敛 性 和 奇 
性 问题 . 在 球 坐 标 下 每 个 球面 区 域 为 


(£-5«e« t ee)n( -isos us). 


其 中 , 6 是 一 个 依赖 于 网 格 步 长 的 小 量 以 确保 6 个 网 
格 之 间 有 重合 . 每 一 个 计算 区 域 均 由 式 (31) 定义 , 但 
是 各 自在 不 同 的 球 坐标 系 下 , 每 一 区 域 在 相应 坐标 系 
空间 是 一 个 矩形 , 6 个 网 格 之 间 的 坐标 可 以 相互 转化 ， 
对 6 个 网 格 区 域 的 详细 描述 可 以 参考 文献 [28]. 
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图 5 AUSM 三 种 算法 vortex 算 例 得 到 在 
y = x 处 的 压强 值 (t = 3) 
Fig.5 Pressure field of three AUSM algorithms to solve 
the vortex problem at y = x (t = 3) 
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将 控制 方程 式 (25) 变换 到 太阳 共 转 坐标 系 下 ， 
考虑 太阳 引力 , 并 添加 动量 和 能 量 源 项 以 描述 太阳 风 
加 速 和 日 冕 加 热 过 程 , 此 时 控制 方程 变 为 


Op 
re pv=N, (31) 
O(pv) B? B-B) 
-a y poarta ———] = 
pe E Sy, 
T 
(32) 
93 HEV.(v- B- B-o) t Ve - 0, (33) 
QE B(v.B)| 
E ee [ees 2:2) 
GM; 
-ø| ; r + Q(Qr)| Bg (34) 
T 
op 2 E ch 
En +V- B zd (35) 
其 中 ， 
1 2 p l5 
E — 5pv viocg ges 
1 
pw-—p-c 3B^ 
7 为 以 日 心 为 坐标 原点 的 位 置 矢量 , -GM [7 为 太阳 


表面 的 重力 加 速度 , 2 = 22/27.2753rad d! 为 太阳 
的 自转 角速度 , 比 热 Y = 1.05. 方程 中 的 Sw 和 Sp 
为 反映 太阳 风 加 速 和 日 园 加 热 过 程 的 动量 和 人 能量 源 
项 , 取 法 与 文献 [28] 相同 . 
采用 Parker?! 的 太阳 风 解 作为 流 场 的 初始 条 
ff. 底面 温度 T = 13 x 105 K, 底面 数 密度 n = 
1.5 x 108cm “， 采 用 如 下 偶 极 磁场 和 8 级 磁场 的 
线性 组 合 
B,—aBO +bB®, Bo -aB +0B®, 
B® ~ 2R3 cos0/r3, BP = R? sin0/r?, 
B(9 = R5[4(5 cos? 0 — 3)/(3r?)] cos 6, 
B? = R*[(5cos? 0 — 1)/rs]sing (a — 3,0 — 1) 
作为 磁场 的 初始 条 件 . 图 6 给 出 了 初始 磁场 的 位 形 . 
模型 计算 30h 后 达到 稳 态 ， 图 7 给 出 了 极 
区 (L = 90) 和 赤道 (L = 2) 处 太阳 风速 度 和 密度 
随 半径 的 变化 , 可 见 太阳 风 速度 随 半径 增 大 , 而 密度 
随 半径 减 小 , 这 与 观测 和 其 他 模型 的 计算 结果 相符 . 
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图 8 给 出 了 磁场 结构 在 子午 面 内 的 投影 , 可 见 在 高 纬 
EKR, 磁场 开放 向 行星 际 空间 ; 而 在 低 纬 区域 , 则 
形成 了 闭合 的 园 环 结构 .高 纬 与 低 纬 磁场 结构 的 差 
异 通 过 太阳 风 与 磁场 的 相互 作用 反映 到 密度 和 速度 


图 6 初始 磁场 在 子午 面 内 的 投影 
Fig.6 Initial magnetic field in the meridian 


plane of projection 
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图 7 稳 态 后 赤道 和 极 区 附近 随 半 径 变 化 的 速度 和 和 密度 
Fig.7 Velocity and density at the equator and polar 


region with radius under steady station 
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结构 上 , 赤道 处 的 太阳 风速 度 比 极 区 高 , 而 密度 比 极 
区 低 , 这 些 结果 与 文献 [1] 的 计算 结果 一 致 图 9 给 
出 了 计算 区 域内 最 大 的 磁场 散 度 误差 随时 间 的 变化 , 
最 大 磁场 散 度 误差 被 控制 在 10 以 下 , 随 着 计算 进 
行 而 逐步 稳定 , 说 明 HDC JATE H AAH fE 
够 有 效 控制 磁场 散 度 误差 
4 Fe 

XE T AUSM, AUSM-- 和 AUSMPW4- 三 种 算 
法 在 使 用 HDC 散 度 消去 方法 后 计算 多 维 MED 问题 
的 性 能 , 结果 表明 AUSM+ 算法 的 性 能 最 好 . 进一步 
使 用 AUSM+ 算法 在 6 块 网 格 构造 下 模拟 了 日 园 结 
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图 8 稳 态 后 磁场 结构 
Fig.8 Steady magnetic field structure 
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图 9 数值 模拟 过 程 中 计算 域内 磁场 散 度 最 大 值 
Fig.9 Maximum of the magnetic field divergence in 


the process of numerical simulation 
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Tg, YTEAS ARR BGCRESELE EE AE SR TTE HAIK 
尺度 结构 , 同时 HDC 方法 能 够 较 好 地 控制 磁场 散 度 
误差 , 具有 进一步 的 应 用 价值 . 
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